Impurity-induced dephasing of Andreev states 
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A study is presented concerning the influence of flicker noise in the junction transparency on coher- 
ent transport in Andreev states. The amount of dephasing is estimated for a microwave-activated 
quantum interferometer. Possibilities of experimentally investigating the coupling between a super- 
conducting quantum point contact and its electromagnetic environment are discussed. 
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I. INTRODUCTION 

How much can nicker noise in the junction trans- 
parency affect coherent transport in Andreev states 
present in a superconducting quantum point contact 
(SQPC)? The assumption of coherent transport in An- 
dreev states is widely used in theoretical work, see e. g. 
Refs. |] and ||. However, in realistic systems interac- 
tion with a dynamical environment will always introduce 
some amount of dephasing, see Refs. ||and [|for a review. 

In the so-called microwave-activated quantum interfer- 
ometers full coherency of Andreev states is assumed and 
a method for Andreev-level spectroscopy is presented. 
The Andreev levels are probed with a microwave field, 
resulting in an interference pattern in the current. If 
dephasing is present this interference pattern will begin 
to deteriorate. This connection between dephasing and 
current makes the microwave-activated quantum inter- 
ferometer a suitable system to study the effect of flicker 
noise in the junction transparency on transport through 
Andreev states. It also provides an excellent opportunity 
to probe the coupling between a SQPC and its electro- 
magnetic environment. 

In the following we study the role of dephasing induced 
by flicker noise in the normal state junction transparency, 
D, of a SQPC in the transport through Andreev states. 
The concrete system which will be considered is the above 
presented interferometer .d Flicker noise can be caused by 
the presence of an impurity atom close to the junction 
which has two states of almost equal energy to choose 
from. When the atom tunnels between its two states the 
junction transparency will fluctuate. Another source of 
flicker noise is the tunneling of an electron between im- 
purity atoms in a doped region. If there are two such 
neighboring defects with available states, a hybrid two- 
level state is formed and the electron can hop between 
the two. This hopping will then add a fluctuation to the 
junction transparency. The amplitude of these fluctu- 
ations depend on the distance between the defects and 
the junction. From now on we will refer to these dynamic 
defects as two-level elementary fluctuators (EFs). 

The microwave-activated quantum interferometer 
(MAQI) is based on a short, weakly biased SQPC which 



is subject to microwave irradiation. Confined to the con- 
tact area there are current carrying Andreev states. The 
corresponding energy levels - Andreev levels - are found 
in pairs within the superconductor energy gap A, one 
below and one above the Fermi level. When a SQPC 
is short (L -C £o where L is the length of the junction 
while £o is the superconductor coherence length), there 
is only one pair of Andreev levels and their positions de- 
pend on the order parameter phase difference, <p, across 
the contact as 



E± = ±E{<f>) = ±AJl - DskTO/2) 



(1) 



Within this pair, the two states carry current in oppo- 
site directions and in equilibrium only the lower state is 
populated. The applied bias, V, through the Joseph- 
son relation <j) = 2eV/h, forces the Andreev levels to 
move adiabatically within the energy gap with a period 



of T p = hir/eV, see Fig. |. 




<t>(t A ) 



<Kt B ) 2ti <|>(t) 



FIG. 1. Time evolution of the Andreev levels within the en- 
ergy gap of a single-mode SQPC. The wavy lines connecting 
points Ai, A2 and Bi, B2 symbolize resonant transitions be- 
tween the levels induced by an applied microwave field. The 
symbol £ defines the position of the resonances. 

The microwave field induces Landau-Zener (LZ) tran- 
sitions between the Andreev levels (symbolized by wavy 
lines in Fig. |l|). If the upper level is populated after the 
second transition a delocalized quasiparticle excitation 
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will be created when this Andreev level merges with the 
continuum. The result will be a dc contribution to the 
current. Further, this current exhibits an interference 
pattern since there are two paths available to the upper 
level. It is this "interference effect" which is utilized in 
the MAQI for Andreev-level spectroscopy. 

In order to estimate the influence of flicker noise on 
the interference pattern we continue with a quantitative 
presentation of the interferometer and the model used for 
the fluctuations. 



II. MODEL AND METHOD 



where r 2 is the LZ transition amplitude, which depends 
on the bias voltage and the amplitude of the perturba- 
tion (0 is the phase of the LZ transition, which can be 
considered constant ).B The phase 4> which inhibits the 
interference is calculated through 

♦4£*(™-t)* < 6 » 

Flicker noise of the junction transparency, D, which is 
the subject to be discussed now, enters through E(4>), 
Eq. (|l|), in this expression. 



Consider a short single-mode SQPC which is subject 
to a high frequency microwave field (hu = 2E((f>) < 2A). 
Let the contact, placed at x = 0, be characterized by 
an energy-independent transparency D. A weak bias, 
eV <C A, is applied across the junction. We choose to 
describe the quasiparticles in the contact region with the 
following wave function, 



V(x,t) =u+(x,t)e lkFX +u-(x,t)e 



-ikpx 



(2) 



where the envelope functions u±(x,t), left and right 
movers, are two-component vectors in electron-hole 
space. To simplify notation we introduce the four- 
component vector u = [u+, u_]. This vector satisfies the 
time-dependent Bogoliubov-de Gennes equation 

ihdu/dt = [Ho + <Jz V g (t)]u , (3) 
Ho = —i~hvF<J z T z d/dx + A [a x cos(0(i)/2) 

+ sgn(xK sin(0(t)/2)] , (4) 

where <7j and Tj denote Pauli matrices in electron- 
hole space and in ± space, respectively, while V g (t) = 

cos(W) is the time-dependent gate potential. We as- 
sume eV^ <C A. 

The boundary condition at x = is 

u(+0) = D- 1 ' 2 \l - t v (1 - D) 1 ' 2 ] u(-0) . 

As a result of the applied high frequency field there 
will be resonant transitions between the Andreev levels. 
These transitions introduce a mixed state which can be 
described within the resonance approximation as 



=pwt/2 



u(x,t) =^2b ± (t) u ± (x)e 



where u + (x) and u~ (x) are the envelope functions of the 
upper and lower Andreev states, while b + and b~ are the 
corresponding probability amplitudes. The final result is 
a dc current through the SQPCjil 



Idc = 



e(2A - Tluj) 



l>-[^)\ =->W 2 (B <!>). 



A. Flicker noise in the junction transparency 

The main topic of this work is to study the effect of 
fluctuations in the junction transparency on the MAQI. 
For simplicity we choose to model the sources of these 
fluctuations, the EFs, with the so-called random tele- 
graph process. This process is characterized by a random 
quantity which has the value +1 or —1 depending 
on whether the upper or lower EF state is occupied. We 
assume that the probability of each state is the same, 
namely 1/2. This is acceptable since EFs with inter- 
level distances, Ei, smaller than k^T, will be "frozen" 
— they behave as static impurities which do not affect 
the dynamic fluctuations of D. In this model the EF 
switches between its two states randomly in time. Physi- 
cally, switching is a result of interactions between the EF 
and phonons or electrons in the contact area. 

In the presence of EFs the junction transparency will 
be modulated. In other words, D — > D + Df(t), where 
Df(t) is assumed to be small. Generally 



with Aj being the coupling strength of the i-th EF. We 
assume that the random processes in different EFs are 
not correlated. Consequently, after a change in variables 
from t to <fi we can specify the random telegraph processes 
(t) through the correlation function, 



,—27»|02- 



(7) 



where 7$ is the switching rate of i-th EF in units of the 
Josephson frequency wj = 2eV/h. It is related to the 
dimensional switching rate T as 7 = hT/2eV. 

Depending on the construction of the SQPC there can 
be any number of EFs which are "in range" to influence 
the transparency. In junctions which are very small it is 
probable that only one single EF will be in the vicinity 
of the contact. In this case, the coupling constant A and 
the switching rate T can be directly evaluated from the 
measured telegraph noise intensity in the normal state, 



I = (2e/hir)r 2 {l-r 2 ){2A-huj) , 



(5) 



S(T) = (I(t + r)I(t)) t 



(I(t))l 
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Indeed, the current through a single mode QPC at low 
temperatures can be expressed according to the Landauer 
formula as I(t) = 2e 2 VD(t)/h. Consequently, the ran- 
dom telegraph noise intensity is equal to 



S(t) = (2e 2 VA/h) 2 exp(-2IY) 



(8) 



and both A and V can be extracted from measured S(t). 
A possible approach for extracting these model parame- 
ters from noise measurements in the case of many fluc- 
tutators in the QPC area will be discussed later. 



III. SINGLE EF 

In the case of a very small contact it is possible to 
consider only one EF and put i = 1. We start by decom- 
posing Eq. (|l|) as 

E{cf>) 

E± = 



E+{4>) + 



1 



-[£(#D0±£(#D-i)], 



(9) 



where £>±i are the two different values the transparency 
fluctuates between. Further, we assume that both 7 and 
gf s = E-/eV are much smaller than the reduced inter- 
level distance (E + ss A)/eV. This means that all devi- 
ations in time are much longer than the Andreev level 
formation time, which is of the order h/2A. 

Fortunately, the expression above, Eq. (^), is linear in 
£ and we can write the effect of the EF as an additive 
contribution to the accumulated phase, Eq. (^]), without 
making any approximations. Namely, $ — * $ + <I> / , with 



9 Mm d<t>. 



(10) 



<t>A 



with gf s (4>) — E-(4>) / eV . After averaging over the real- 
izations of the random process the expression (|^) 
for the MAQI current is replaced by 



he =h[l-W cos (29 + 2$)] 



(11) 



where W = (e f ) contains the dephasing. Without 
dephasing W = 1 and the effect of the phase &(4>a, 4>b) 
is a modulation of the dc current between and 2 Jo, the 
modulation depth 



(h 



/mm) = \W\ 



(12) 



being equal to one. This modulation of the current is the 
interference effect utilized in the MAQI. When dephas- 
ing enters, the quantity W(4>a, 4>b) will decrease and the 
modulation \W\ of the current envelope, i. e. the inter- 
ference pattern, will in turn decrease. 

To facilitate the calculation of the dephasing term, W, 
we define the auxiliary function, 



*(0) 



exp 



i / W gfM'W) 



4>A 



(13) 



The quantity of interest, W, is related to as W = 

^((ps)- Further, it can be shown, cf. with Ref. ^, that 
Eq. (O) satisfies the differential equation 



d 2 * 



2 7 - 



dlng /s (0) 



— + ff/,fa)* = 0, (14) 



with the initial conditions 

*(^a) = 1, d*/d^| 0=0A =0. 



(15) 



Let us consider the following two cases in more detail: 
(i) the "slow EF", 7 -C <?/ s , which corresponds to low 
temperatures, and (ii) the "fast EF", 7 ^> gj s , which 
corresponds to relatively high temperatures. 



A. Slow switching EF 

In the low-temperature limit the EF will slowly switch 
between its two states, and we can let 7 — * in Eq. (pH). 
Further, by introducing the function, 



S(0) = / tygfsMMt), 



(16) 



and applying the initial conditions (|l^) we obtain the 
integral equation, 

S 2 + * 2 = 1 . 



The solution follows as *^(^) = cos $(</>) where $(0) 
jf d(j>' gf s (4>')Y Consequently, = cos$^ with 



r<t>B 

$ c = = / d<t>g fs (<t>) 

J <p A 



(17) 



A better approximation can be found by looking for the 
solution in the form ^>^((p) = u((p)^/^' ((f>) and assum- 
ing u{4>) to be a slow function. Neglecting the second 
derivative of u((f>) we obtain the differential equation 



du/d<t> — —7 r)(<f)) u 




dg fs 



Defining 



2 ff 2 s tan$(0) 



^ r<t>B 

2 S> J<t>A 



(18) 



(19) 



(20) 



we arrive at the following expression for the oscillating 
part of the current, 

W cos($) = e - 2lv ^ cos($ c ) cos($) 

= e - 2 ^ c [cos($ + $ c ) + cos($ - $ c )] /2 . (21) 
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For a constant gf s , vq = \ and = 2g/ s £, where £ = 
(0b — 0a) /2 is equal to half the distance between the 
resonance positions, see Fig. |l|. In the general case these 
quantities are increasing functions of £. 

At 7 = 0, the current is split into two interference 
patterns of equal magnitude shifted by the phase $f. A 
plot of as a function of £ is presented in Fig. g. The 
transparency D and the strength of the fluctuator are 
shown in the inset. There is no dephasing, only a dis- 
tortion of the interference pattern, the modulation being 
|W| = |cos$ f |. 

This splitting into two patterns of equal magnitude 
follows from the assumption that the occupation prob- 
ability is the same for the two EF states. The general 
case of arbitrary probabilities for the EF states can be 
solved numerically. At finite 7 dephasing takes place and 
the amplitude of the interference oscillations decreases by 
e~ 2lv ^ . The physical reason of dephasing is the finite life 
time of an EF in a given state. We have calculated W for 
7 = 0.1, comparing the analytical approximation above 
against a numerical solution of Eq. (14), cf. with Fig. [|. 
(The plot for the analytical case when D — 0.6, A = 0.05 
is missing because of numerical difficulties) . 



D=0.7, A=0.1 

- - ■ D=0.6, A=0.1 

— - D=0.6, A=0.05 




FIG. 2. In the case of flicker noise from a single EF in 
the low temperature limit the MAQI interference pattern is 
split into two parts, phase-shifted by shown here as a 
function of the resonance position f. The inset shows the 
contact transparency and the strength of the fluctuator. Here 
eV/A » 0.1. 
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FIG. 3. In the case of flicker noise from a single EF at low 
temperatures, when the hopping rate 7 is non-zero, there will 
always be some amount of dephasing (see text). This figure 
shows the dephasing factor W as a function of the resonance 
position £ when 7 = 0.1. The two curves marked with dia- 
monds (o) are calculated from W in Eq. (pl|), while the rest 
are direct numerical solutions of Eq. (141). The inset shows 
the contact transparency and the strength of the fluctuator. 
In this plot eV/A « 0.1. 



B. Fast switching EF 

In the case of fast switching the differential equation 
( |l4| ) can be approximated as, 

.,2 



d* 

d0 



g}M 
27 



The solution in this case is easily obtained, as 
W = = exp[-K(<t>BM], 

with 

1 



K(cf> B ,cl> A ) = — dcjjgfM- 
z 7 J 4> A 



(22) 
(23) 
(24) 
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- ■ D=0.6, A=0.1 
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FIG. 4. In the high temperature hmit, flicker noise of a 
single EF leads to damping of the MAQI interference pattern 
by a factor of W — exp(-K). The inset shows the contact 
transparency and the strength of the fluctuator. In this plot 
eV/A ~ 0.1 and the switching rate is V = 3. 

Here we find an exponential decay of the interference 
term. Expression (24) describes an effect which is sim- 
ilar to motional narrowing of spectral lines .0 When the 
EF fluctuates rapidly enough comparied to the "energy 
resolution" E_/h, influence from the difference between 
two EF states is smeared and dephasing will be of a dif- 
fusive character, with en effective, time-dependent, dif- 
fusion constant g'j s /2^. A calculation of the factor K for 
7 = 3 is shown in Fig. |] for three different sets of val- 
ues of the junction transparency and the EF's strength. 
The effect of dephasing from flicker noise will decrease 
for higher bias voltages, since a higher bias, through the 
Josephson relation, leads to a shorter time between the 
resonances. 

To get an idea about the general case of arbi- 
trary switching rates, it is instructive to consider a (f>- 
independent gj s . This approximation is only valid when 
\dE ((/))/ d(j>\ is small, but should provide the general be- 
havior of W. In this limit Eq. (Q) can be solved analyt- 
ically, resulting in 



W 



-2 7 C 



r 



9% 



:>i"l:(2 M J 2 9% 



(25) 



From this expression one can conclude that there will al- 
ways enter an exponential decay of the interference cur- 
rent, except for the limit of 7 — > 0. In this limit there 
should only be an additional oscillation added to the cur- 
rent as a function of £ since for 7 < |<7/ s |, the hyperbolic 
functions behave as trigonometric functions. 



IV. MANY EF'S 

Let us consider a large number of EFs with varying 
switching rates distributed in the contact area. For sim- 
plicity, we shall assume that only the fluctuators with 
inter- level spacings Ui < kT are important, and that 
their distribution is uniform, Vu{U) = PqV. Here V is 
the sample volume. Further, we assume that the switch- 
ing rates 7$ are the same for both transition directions 
(up och down) between the EF's levels. This assumption 
is natural because the ratio between the corresponding 
transition rates is exjp(—Ui/kT). Within the assumptions 
discussed above, the final results are substantially simpli- 
fied while preserving the essential depedence on temper- 
ature and the resonance position. These approximations 
agree with a general theory developed in Ref. || for the 
case of dephasing by two-level systems (TLS) in glasses. 



The first step now is to linearize the SQPC's trans- 
parency with respect to & as D — > D + Df(t). This 
allows us to once again find an additive contribution to 
the accumulated phase, Eq. (|^), which in this case will 
be, 



4> 



5> 



4>A 



&(</>) 9fm{4>) 



(26) 



Here we have defined gfm{4>) = {l/2eV)dE{4>)/dD. In 
the same manner as in Sect. QT0, but this time averaging 
over A's, 7's and £'s, we can express the modified MAQI 
current through expression (xTji) with 



W 



&(0)d0 



(27) 



To approximate this average we use the Holtsmark 
methocf] which is valid in the limit of many fluctuators, 
N = P V kT > 1. This allows us to rewrite Eq. @ as 
the average over the contributions from single EFs, 



W s (A,7) = (exp [iA g fm (4>)m 




as 



W « exp (-PoVfcT (1 - W S (A, 7)) 



A, 7 



(28) 



(29) 



Since the number of EFs is assumed to be large, to keep 
dephasing at a reasonable level it is important to keep 
(1-Ws) small. 

With the solutions for W s found in Sect. Ill the av- 
erage (1 — W s ) remains to be calculated. To calculate 
this average one has to specify the distributions of the 
parameters A and 7. The simplest and most natural 
assumption is that these two quantities are not corre- 
lated. Consequently, the distribution V(A, , y) can be de- 
coupled as Va{A) V 1 ( r y). To specify the distribution Va 
let us assume that the EFs are uniformly distributed in 
space. An EF behaves like a dipole, either electric or elas- 
tic, this allows us to specify it's interaction strength as 
A(r) = A /r 3 , where^r is the distance between the con- 
tact and a given EF,EHl while Aq is a coupling constant de- 
pendent on a specific interaction mechanism. Note that 
the quantity Aq has dimension of volume. Within this 
model we arrive at the normalized distribution function 
V A (A) = AttA /3VA 2 (see Appendix B for details). The 
distribution 7^(7) is specified in a manner which is com- 
monly used in glasses. Namely, the logarithm of 7 is as- 
sumed to be uniformly distributed. Hence, ^(7) oc 7 _1 , 
see Appendix [b|. To normalize it let us take into account 
that for a given energy spacing U there is a maximal 
switching rate. Since we are interested in the fluctua- 
tors with Ui < kT, we can specify the maximal switching 
rate as jt, which is a function of the temperature. The 
actual temperature dependence is determined by the spe- 
cific interaction mechanism between the EF and its en- 
vironment. If the transitions between the EF states are 
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caused by interaction with phonons, then oc T 3 ,0 
while if the transitions a*e caused by the electrons ex- 
citations, then 7x oc T.H3 Therefore , the normalized 
distribution can be specified as "P 7 (7) = (£7) , where 
L = ln(7T/7mm) 1- Here we have introduced the min- 
imal switching rate, 7 m i n - To express the decay in a more 
clear form let us introduce the dimensionless frequency 
VcL corresponding to the interaction strength for an EF 
separated from the contact by an average distance to the 
active fluctuators, f = (AirPokT/S)^ 1 / 3 , divided by the 
Josephson energy 2eV. We can specify v d as 



v d = 4irP kTA /3 = A /f 3 



(30) 



The decay rate /C 
sion 



In W is then given by the expres- 



K. 



C 



dA 
A* 



, 7 



[l-W a (A, 7 )] 



(31) 



Note that gj s has to be replaced with Agf m in the expres- 
sions for W s in the many EF case, because of differences 
in notation. 

To estimate the amount of dephasing let us take into 
account that the asymptotic expressions for W s , found 
in (HJ) and ( ^3|) for slow and fast switching respectively, 
match at 



7U C w (Ag ( ) 2 r} C /2-y , 



(32) 



where vq differs from defined in Eq. ( pp| ) by the re 
placement gf s — > Agf rn in the differential equation (18) 
while 



9( = / !lf mi") • </,, 



1 



C J <t>A 



g) m {^)M- (33) 



Note that vq is a function of A, v^(A). Consequently, 
Eq. ( |32[ ) should be treated as an equation to determine 
the characteristic value of the coupling A. Defining the 
solution of Eq. (j32| ) as Aq and splitting the integration 
over A in Eq. (|31[) as 



W.(A,j)\ (34) 



we can use the expression fl23| ) in the first interval and the 
expression (^lj) in the second one. Since both integrals 
are determined by A^ we arrive at the estimate 




-InWw3i/ d C/c. 



(35) 



where = g^y/r]^v^(A^) ^ s some function of £, rather 
smooth if D is not close to 1. We see that the interference 
pattern decays exponentially with an increasing distance, 
2£, between the resonances. Generally the ^-dependence 
of can be calculated numerically for a given trans- 
parency D using the analytical expressions (20) and (32). 
For a constant gf m , fa = gfm- We do not analyze here 
the function in detail. 



A. Non-optimal EFs 

In the previous sections it has been assumed that the 
system size is infinite. A consequence of this assumption 
is that, independent of temperature, the EFs which have 
the strongest effect on the junction transparency will al- 
ways be included in the estimates. From the method used 
to obtain the general estimate for many EFs, Eq. (|3^), 
one can conclude that the EFs which have the most ef- 
fective coupling fulfill A this corresponds to a spa- 
tial distance r 7 = (^4 Mc) 1/3 ~ (^o5c/t) 1/3 - A further 
point is that the rate 7 is confined to the interval between 
7 m i n and "fx- Thus we have actually assumed that the 
size of the region where EFs reside has a size greater than 
r max w (Ao i g^/7 m in) 1 ^ 3 , and that there is no "excluded 
region" without EFs near the contact with the size less 
that r min w (Ao<? c /7t) 1/3 . 

What happens if this "optimum" EF is out of the 
range? This can occur if the system is limited in size, 
or if there is a specifically pure region around the con- 
tact. 



1. Role of finite size of the sample 

Let us first discuss the role of a finite size, R , of the 
region containing EFs. If R < r m i n = (Aq/A^) 1 / 3 then 
all EFs will act as "slow" ones. To estimate In IF in 
this case one can use the expression W s = cos for the 
whole integration region over A in Eq. ©, arriving at 



\nW = 2is d g c (F(2Aog c (/R 3 ). 



Here 



F(z) 



-dx . 



(36) 



(37) 



which is a decreasing function of its argument. This ex- 
pression (|3^) is valid if its right-hand side is less the right- 
hand side of Eq. (§|). 

At intermediate values of R, 

(A)ff C /7min) 1/3 » R » (^o.9 C /7t) 1/3 , 

as in the previous case, only the second integral in the 
expression (|34|) does exist, lowest limit should also be 
replaced by A^g^/ R 3 . However, the approximation W s = 
cos is not valid any more. Using the approximation 
(pl|) one can obtain 



W = exp 



- v dC,fc 



\n(-f T R 3 /A g c ) 
ln(7r/7min) 



(38) 



This means that when the size of system is limited the 
amount of dephasing can be less than estimated for an 
infinite system. 
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2. Role of the spacer 

Let us now discuss the role of a "pure" region (spacer) 
near the SQPC where there are no EFs. If a typical 
diameter ro of such region is large enough, such that 

1"0 ^ fmax = (^Offc/Tmin) 1 / 3 , tnen aU EFs act as "fast" 

ones. 

The single EF solut ion, W s , in this limit is given by 



Eq. 

Eq. (29j) we arrive at 



in Sect. IIIB. After calculating the average in 



InW 



(39) 



It is difficult to estimate the actual amount of dephasing. 
We have to restrict our conclusions to interpreting how 
the amount of dephasing will change depending on the 
parameters 7t, 7min, C an d T. We already know that 
dephasing will increase with temperature. Generally one 
can also say that dephasing will increase with £. How- 
ever, at large enough temperatures when 7 m ; n appears 
large enough, the dephasing will decrease. The free pa- 
rameter of the theory, Aq, can be estimated only roughly 
through comparison with the noise measurements in the 
normal state. To map the parameter Aq to the poise 
let us employ the theory of flicker noise in a QPQL3 to 
the case of a single mode contact. According to that 
theory, results for the noise intensity S(t) are substan- 
tially dependent on the relationship between the maximal 
and minimal distances between the EFs and the QPC. 
The simplest case, which is quite realistic, is when these 
distances are of the same order of magnitude. When 
T^ 1 <C |t| <C TjT^ the noise intensity can be expressed 
as, cf. with Rcf. jl3| , 



S(r) 



2e 2 V\ 2 fATrP kTA \ 2 ln(l/r D 



"I) 



in(r T /r D 



(40) 



By obtaining estimates for / m j n from noise spectra in 
the normal state one can, in principle, estimate the cou- 
pling parameter Aq. A key point is to make measure- 
ments of both the MAQI interference pattern and the 
normal-state noise spectra in a rather large frequency 
range. This combination does not look too simple. 



V. CONCLUSIONS 

We have presented a method for investigating the in- 
fluence of flicker noise in the junction transparency of a 
SQPC on coherent Andreev states. This is done by es- 
timating the effect of these fluctuations on the so-called 
microwave-activated quantum interferometer (MAQI)tl 

For a small contact when only a single EF is in range to 
affect the junction transparency there can be either a dis- 
tortion or a decay of the MAQI interference pattern. A 
distortion appears for very slow switching rates of the EF 
and a weak decay (dephasing) for fast rates. It is possible 



to confirm our model in the fast switching limit experi- 
mentally. The only unknown parameter is the switching 
rate 7 which can be found by measuring the telegraph 
noise of the contact in the normal state. This is best 
done by driving the system into the normal state with a 
magnetic field, since 7 is temperature dependent. When 
7 is known it is then possible to calculate the amount of 
dephasing and compare with experimental results. 

In the limit when the influence of many EFs has to be 
considered we have arrived at more general results. It 
is not possible to make any exact predictions since the 
distribution and coupling strength of the EFs are sample 
dependent. However, our calculations show that in the 
presence of many EFs there will always be an exponential 
decay of the MAQI interference pattern. The strength of 
this dephasing should be about the same for all switching 
rates. One exception is when there is an impurity-free re- 
gion near the SQPC, in this case dephasing will decrease 
for higher rates. 

Finally, we note that this paper together with work 
in Ref. |l4j presents a framework which can be used to 
investigate the coupling of a SQPC to its electromagnetic 
environment. 
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APPENDIX A: DERIVATION OF EQUATION 
FOR W 

In this appendix the differential equation, Eq. [l4[ for 
W is derived. The following is valid for the random tele- 
graph process £(£), 



(Al) 



when 0i > <j>2 > • • • > 4>n- Further, all averages with an 
odd number of £:s are equal to zero. This follows from 
the assumption of equal probability for both EF states. 

To facilitate the calculation of W, the auxiliary func- 
tion 



exp 



4>A 



(A2) 



is defined. By expanding the auxiliary fu nctio n in a Tay- 
lor series and utilizing the relation in Eq. Al above, it is 
possible to rewrite the auxiliary function as 



7 



d0! / dfcgfMgfsi^e- 2 ^^^ 



A derivation with respect to 
integro-diffcrcntial equation 



After normalization we have, 
(A3) P( 7 ) 



d* 

~~7T = -9fs{4>) 



leads to the following 
' ff/s (0')e- 2 ^'-^(<A (A4) 



<j>A 



and a second derivation provides the following differential 
equation 



d 2 * 



27- 



dln gf a {<j>) 



with the initial conditions 

tf(^0 = l, d¥/<ty|^ A =0. 



(A5) 



(A6) 



This equation makes it possible to find approximate an- 
alytical solutions for W in the single EF case. 



APPENDIX B: DISTRIBUTIONS, 7> 7 ( 7 ) AND 
Va{A) 

Here we outline the derivation for the distribution 
functions V (7) and V(A) which are necessary when av- 
eraging in the many EF case. Let us begin with V(A), 
where A is the coupling strength between an EF and 
the contact. By assuming that an EF can be modeled 
with a dipole field, we can state-that A cx An/r 3 where 
r is the distance from the EFE!J The normalized dis- 
tribution is then found through the following integral, 



V(A) = V' 1 JdrS{A~ A /r 3 

P(A) = (4tt/3V) (A /A 2 



(Bl) 



To find the distribution of 7 a closer look at the struc- 
ture of the two level system (TLS) which is assumed for 
the EFs is necessary. The Hamiltonian for a TLS can be 
H = Aa z —Aa x , where A is the energy difference between 
the minima of the two states and A is the tunneling cou- 
pling between the two states. After diagonalization the 
excitation energy is found to be U = -\/A 2 + A 2 . By as- 
suming that the tunneling coupling decays exponentially 
we have that V(A) oc 1/A. Further, the hopping rate 7 
depends on (A/U) 2 and a term U 3 /hE 2 where the latter 
term comes from assuming that hopping is phonon me- 
diated (E c is the parameter characterizing the coupling 
energy between the phonons and the EFs). Within these 
assumptions V(j) is given by, 



^7(7) 



dAdAP{A) 5{U 2 
A 2 U 3 



A 2 ) 



xS 7 



(B2) 



1 



(B3) 

being limiting values 



where C = ln(7x/7roui) with 7^, 
for the EF switching rate 7. 

In the general case V = V(U, 7,^4) and the number 
of EFs will be given by N = J dU djdAV{U, 7, A). By 
assuming a constant distribution for U, which we label 
To, we arrive at the number of EFs as N = PoVkT. 

The final distribution function for A and 7 is then 



V{A, 1 ) = 



AirPnkT An 1 



A 2 £7' 



(B4) 
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